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ABSTRACT 

The currently favored method for estimating radii and other parameters of transiting-planet 
host stars is to match theoretical models to observations of the stellar mean density p*, the ef- 
fective temperature Teff, and the composition parameter [Z]. This explicitly model-dependent 
approach is based on readily-available observations, and results in small formal errors. Its per- 
formance will be central to the reliability of results from ground-based transit surveys such as 
TrES, HAT, and SuperWASP, as well as to the spaceborne missions MOST, CoRoT, and Kepler. 
Here I use two calibration samples of stars (eclipsing binaries and stars for which asteroseismic 
analyses are available) having well-determined masses and radii to estimate the accuracy and 
systematic errors inherent in the method. When matching to the Yonsei-Yale stellar evo- 
lution models, I find the most important systematic error results from selection bias favoring 
rapidly-rotating (hence probably magnetically active) stars among the eclipsing binary sample. 
If unaccounted for, this bias leads to a mass-dependent underestimate of stellar radii by as much 
as 4% for stars of 0.4 Mq, decreasing to zero for masses above about 1.4 Mq. Relative errors in 
estimated stellar masses are 3 times larger than those in radii. The asteroseismic sample suggests 
(albeit with significant uncertainty) that systematic errors are small for slowly-rotating, inactive 
stars. Systematic errors arising from failings of the Yonsei-Yale models of inactive stars probably 
exist, but are difficult to assess because of the small number of well-characterized comparison 
stars having low mass and slow rotation. Poor information about [Z] is an important source of 
random error, and may be a minor source of systematic error as well. With suitable corrections 
for rotation, it is likely that systematic errors in the method can be comparable to or smaller 
than the random errors, yielding radii that are accurate to about 2% for most stars. 
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The radii of transiting extrasolar planets are 
measured in units of the radius of the host star. 
Thus, to understand the key properties of such 
planets, as well as for a sense of their astrophysical 
context, the masses, radii, and compositions of the 
host stars are of great interest. With the advent 
of efficient searches for transit ing planets such as 



Kepler dWalker et all l2003t iBaglin et all 12007 



Koch et al.l l2004f) . the need for accurate stel- 
lar par ameter esti mates has become acute. Re- 
centlv. ISozzetti et al. (2007L_ B.ako s et al. (2007 



Winn et all (l2007l) Icharbonneau et all (|2007l ). 



Torres et al. (12008) and others have built on 



TrES, HAT and SuperWASP (lAlonso et al.l l2004l: 
Bakos et al.ll2002i 12004 IPollacco et al.l l2006t). and 



especially spaceborne analogs MOST, CoRoT, and 



work bv iMandel fc Agol (|2002l) to develop the 
currently-favored approach to stellar parameter 
estimation, which for the purposes of this paper 
I will term "the method" . This paper aims to 
understand the precision and accuracy that one 
may expect from this method, and, by comparing 
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against other ways of measuring stellar parame- 
ters, to identify the sources and likely magnitudes 
of systematic errors that may affect the method's 
results. 

The method uses the transit light curve to 
estimate the 3 parameters Rp/R^, R*/a, and the 
transit impact parameter b, which measures the 
minimum projected distance between the centers 
of the star and planet in units of the stellar radius. 
Here Rp is the planetary radius, i?* is the stellar 
radius, and a is the orbital semimajor axis. The 
fitting procedure used to make these estimates re- 
lies on a model of the stellar atmosphere to char- 
acterize the stellar limb darkening. From i?*/a, 
the planet's orbital period P (which ordinarily is 
known with extremely high accuracy), and Ke- 
pler's third law, one may compute the stellar mean 
density with an accuracy that is limited by the 
photo metric precision of the light cur ve measure- 
ment (jSeager fc Mallen-Orneia3l2003l ). This cal- 
culation depends only on geometry and on Newto- 
nian gravitation, so it is independent of modeling 
assumptions. Next, by comparing observed high- 
resolution spectra of the parent star with model 
stellar atmospheres, in typical cases one can es- 
timate Teff with precision of about 1%, and [Z], 
the logarithmic metal abundance relative to so- 
lar, with quoted errors of about 0.05 dex (e.g. 



Valenti fc Fischer (Hods')). Finally, one searches a 
grid of stellar evolution models for the best match, 
in a sense, between computed and observed 
values of p,, T^s, and [Z]. The search for opti- 
mum model parameters {mass, age, composition} 
is usually conducted using Markov chain Monte 
Carlo (MCMC) methods. These also yield esti- 
mates of the formal distribution of error in (and if 
desired, the covariance between) the various model 
parameters, and indeed any other global property 
of the stellar model (_R», p,, log{g), Tcs, etc.). 

An important virtue of the p, method is that it 
typically yields significantly smaller formal errors 
than do older techniques - uncertainties of 2% or 
less in i?, are commonly achieved. This is possible 
because the observables {p,, Teff, [Z]} are indi- 
vidually fairly well-determined, and because each 
of them corresponds fairly closely to one of the 
model parameters {M,, A,, [Z]} describing the 
stellar mass, age, and metallicity. Thus, during a 
star's main-sequence lifetime. Toff depends mostly 
on Af, , p, depends mostly on A* , and (no surprise 



here) [Z] depends mostly on [Z]. On the other 
hand, the method is manifestly model-dependent. 
To give correct results, the technique depends 
upon having evolution models that accurately re- 
produce the mass-luminosity-radius-composition 
relations followed by real stars. Moreover, the 
method depends upon stellar atmospheres models 
both to account for limb darkening and to allow 
interpretation of the stellar spectrum in terms of 
Teff, log{g), and [Z]. There are good reasons to be- 
lieve that stellar evolution and atmosphere mod- 
els in fact represent real stars fairly well, but one 
nevertheless desires direct test s of the p<, method . 
Fortunately the recent update ( Torres et al. 2009f ) 
of Andersen's (1991) classic review paper on fun- 
damental parameters of stars in eclipsing binary 
(EB) systems collects a large sample of stars with 
well-determined masses and radii. There is also an 
almost independent but far less numerous group of 
stars for which asteroseismology (often combined 
with other measurements) provides similarly pre- 
cise data. Both of these samples of stars are well 
enough characterized to permit tests of the behav- 
ior of the p* method. 

Torres et all (|2009t ) (henceforth TAG) contains 



the results of analyses of 94 EB systems plus a 
Cen A and B, comprising 190 individual stars for 
which masses and radii are thought to be known 
with accuracies of better than 3%. The masses 
and radii thus determined are almost independent 
of complex theory, the only exception being the 
weak dependence of inferred radii on assumptions 
about the stellar limb darkening. Also given by 
TAG are estimates of Toff and, for 21 of the sys- 
tems, [Z]. For most of these systems Toff and [Z] 
were determined from photometry, using meth- 
ods that are not consistent across the sample. 
The spectroscopically-derived estimates of these 
parameters that are common in the exoplanet con- 
text are difficult to obtain for EBs, because of the 
blending, broadening, and relative Doppler shifts 
that characterize double-line spectra. 

From the above information it is straightfor- 
ward to synthesize the input data (p*. Toff, and 
if only by assumption, [Z]) that would be mea- 
sured if a transiting planet were to orbit any of 
the stars. One may then apply the p* method to 
each EB component, and compare the masses and 
radii that emerge from the model-fitting process 
to those that were actually measured. 
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The stars with asteroseismic data consist al- 
most entirely of ones with roughly solar mass, 
and largely of post-main-sequence objects having 
greater than solar luminosity. The observational 
demands for asteroseismology are severe, so in 
most cases the only reliable observable is the so- 
called "large frequency separ ation" At^, which it - 
self depends mostly upon ( Hansen et al. 2004 ). 
In order to obtain a well-constrained result, the 
measured oscillation frequencies are often aug- 
mented with other kinds of data, such as interfer- 
ometric estimates of radii, or dynamical mass es- 
timates in the case of binary systems. Finally, de- 
tailed interpretation of this information always in- 
volves comparisons between stellar evolution and 
oscillation models and the observations. Thus, 
the asteroseismically-determined stellar parame- 
ters are likely no more accurate and are certainly 
as model-dependent as those obtained from the 
method. They usually depend on different model 
implementations, however, so at worst they tell us 
something about consistency among extant mod- 
els of stellar evolution. This in turn provides some 
insight into the poorly-understood possibilities for 
systematic error in the asteroseismic analysis. 

I implemented the method, using the Yonsei- 
Yale evolution tracks (lYi et al. _ 2001 Kim et al 



2OO2I : lYi et al.ll2003l : [Pemarque et al.ll2004h as the 

needed stellar evolution models, and I then applied 
it to the TAG tabulation of EBs and to 15 stars 
with asteroseismic measurements. In the rest of 
this paper, I present the computational methods I 
used and the results I obtained. Section 2 briefly 
describes methodology and algorithms. Section 3 
gives the results of the comparison, illustrated in 
various ways. I find that the p* method applied 
to EB components generates errors in both mass 
and radius that are small but significant, and that 
depend systematically on the basic stellar param- 
eters. Finally, in Section 4 I investigate the origin 
of these systematic errors. 

2. Parameter Estimation Methodology 

Estimating stellar parameters first requires the 
observations that are to be fitted, and their un- 
certainties. I took the observables to be p*, Tcs, 
and [Z]. I derived uncertainties for the latter two 
quantities in the obvious way, from the stated un- 
certainties in the photometric data. The uncer- 



tainty in p* required a different treatment, how- 
ever. My chief purpose is to assess the system- 
atic errors committed by the p* method because 
of failings in modeling stellar structure. For this 
purpose it makes sense to treat the EB and aster- 
oseismic masses and radii as error-free, and then 
assign an uncertainty to p, that is small enough 
that its effect on the derived parameters is unim- 
portant. Thus, I took the uncertainty in p* to be 
about 4.5%, which is typical of what might be ex- 
pected from a very good ground-based observation 
based on several transits, or a rather poor space- 
borne one. 

The central feature of the p* method is the 
stellar evolution model. I used the Yonsei-Yale 
(henceforth YY) evolution tracks downloaded 
from their "Evolutionary Tracks" site. These 
models span masses in the range 0.4 M© to 5 M©, 
and metallicity between 0.001 and 0.08 {[Z] be- 
tween -1.230 and -1-0.673, assuming a solar metal- 
licity of 0.017). They also employ an initial he- 
lium abundance Y — 0.23 + 2Z, and a constant 
mixing length of 1.7432 times the pressure scale 
height. The temporal evolution of each model is 
tabulated on a nonuniform grid of stellar age A*, 
with the same number of steps for each stellar 
mass, and a given step number corresponding to 
the same evolutionary state (meaning core helium 
abundance or core mass) for every stellar mass. 
This grid based on evolutionary state is conve- 
nient for interpolation between different stellar 
masses, since one need not do a search in the age 
dimension to locate models of different masses 
that have otherwise similar properties. The YY 
model grid includes models for which the abun- 
dances of high-alpha elements (O, Ne, Mg, etc.) 
are enhanced relative to the solar composition. In 
this study I have not investigated the influence of 
this degree of freedom, however. 

To apply the YY model grid, I first used the 
interpolation routines supplied with the models to 
resample them onto a uniform grid in logAf, and 
[Z]] for the computations described here I used 
150 mass steps between 0.4 Mq and 5 M©, and 25 
steps in [Z] between -1.230 and 0.673. I did not re- 
sample the tables in the age dimension. From the 
resamplcd grids giving Af», i?,, and the luminos- 
ity i», I then precomputed similar tables giving 



-'http: / /www. astro. yale.edu/domarque/yystar. html 
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other quantities of interest, such as and T^ff. 
Finally, to obtain parameters at arbitrary values 
of {A*, M*, [Z]}, I interpolated into these resam- 
pled grids, using cubic interpolation in the age- 
mass plane and linear interpolation in [Z]. 

The MCMC algor ithm has been described by 
Tegmark et al. ( 2004 ) , a nd in the con t ext of tran- 



siting planets, by e.g 



Charbonneau et alj (|2007r ). I perform the Markov 



Winn et all (|2007[ ) and 



chain random walk in the space of indices into the 
3-dimensional grid of stellar evolution models. I 
use a Metropolis-Hastings algorithm with Gibbs 
sampling and random permutations of the order 
in which the indices for age, mass, and [Z] are 
varied. The probability distribution for accept- 
ing a step that decreases the merit function was 
Gaussian, with standard deviations along the 3 
axes chosen so that steps along each axis succeed 
about 25% of the time. In addition to estimates 
of the mean values and standard deviations of the 
model parameters, the software produced many 
diagnostics. These included marginal probability 
distributions for every parameter of interest, scat- 
ter plots, and various convergence measures, so 
that misbehavior of the MCMC algorithm or am- 
biguities in the models could be identified. I found 
the most useful convergence test was visual inspec- 
tion of the parameter chains, the corresponding 
chain of values, and various 2-dimensional scat- 
ter plots of these quantities. In doubtful cases I 
reran the MCMC process using chains that were 
lengthened by factors of 2 and 4; these longer runs 
always yielded chains that were consistent with 
each other, and that gave convincingly smooth 
and well-sampled sample distributions in parame- 
ter space. 

The merit function itself was intended to rep- 
resent the posterior probability attached to each 
set of stellar model parameters, taking into ac- 
count the observations. From a Bayesian perspec- 
tive, the logarithm of this probability is the sum 
of — x^/2, representing the probability that the 
model is consistent with the observations, and the 
logarithm of a prior probability, representing the 
distribution of model probabilities in the absence 
of any observational data. Each of the parameters 
{A^,, M*, [Z]} has a prior probability associated 
with it, but these are justified in different ways for 
the various parameters. 

An implication of the time grid chosen for the 



models was that the interval (in years) correspond- 
ing to a single time step was highly variable, both 
for a single choice of mass and composition, and 
from one such choice to another. The variations 
in this probability are artifacts of the model im- 
plementation and can be quite large, so it seemed 
important to account for them. Therefore I took 
the prior probability for each model at each age to 
be proportional to the duration of the correspond- 
ing time step, normalized by the maximum age 
of the star represented by the model. Results of 
the MCMC process seldom depended strongly on 
whether I imposed the age-dependent prior prob- 
ability. Exceptions to this rule occurred for stars 
having similar radii and temperatures at two dis- 
tinct ages, one before and one after arrival on the 
main sequence. In these cases the age prior proba- 
bility tended to discriminate against the pre-main- 
sequence phase, because it is very short-lived. 

One could also impose prior probabilities based 
on the statistical distribution of stars (by mass 
and by composition) in the solar neighborhood. I 
elected not to do this, since the detailed choices 
of priors in these cases are somewhat arbitrary, 
and in any case the observed values of mass and 
[Z] are so tightly constrained that reasonable prior 
probabilities have little effect on the calculation's 
results. 

Although the MCMC procedure gave a robust 
estimate of mean values and corresponding dis- 
persions, it was an inefficient way to determine 
the exact position of the merit function optimum. 
But this information was valuable for estimating 
sensitivities (eg to assumed biases in Tcff or [Z]), 
which must be calculated by measuring the change 
in optimum values for small changes in the obser- 
vations. For this reason I also computed the pa- 
rameters giving the optimum merit function using 
an "amoeba" search, which when started near the 
best parameters found by MCMC, converged to 
the local optimum with negligible error in a few 
tens of iterations. In what follows, unless other- 
wise stated, I use the amoeba-search values for es- 
timated parameters, with uncertainties computed 
from the MCMC results. 

The probability distributions generated by the 
MCMC method were usually close to multivari- 
ate Gaussian, reflecting the assumed distribution 
of errors in the input data. But in a substan- 
tial minority of cases, the stellar models permitted 
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multiple solutions that fit the data almost equally 
well. In such cases I encountered the usual prob- 
lems associated with complicated merit-function 
behavior, among them that mean values of param- 
eters found by MCMC did not agree with the opti- 
mum values found by the amoeba search, and the 
amoeba process itself was likely to yield different 
results depending upon its starting values. These 
problems add scatter to the results (particularly 
when computing sensitivities of stellar parameters 
to changing input data), but these uncertainties 
are not large enough to affect the conclusions. 

3. Application to Eclipsing Binaries 

The list of stars supplied by TAG contained 95 
binaries, or 190 stars for which masses, radii, and 
uncertainty estimates were provided. Of these, 34 
had M* > 5M0 or M < OAMq, which are the 
limits of the mass range covered by the YY mod- 
els. This left 156 stars that could be analyzed by 
the procedure. Only 32 of these had individ- 
ual estimates of [Fe/H], which I simply equated 
to [Z]. For a few stars, the estimated [Z] values 
were derived not from direct observation of the 
star in question, but from the star's presumed par- 
ent population (star cluster or galaxy). For stars 
without explicit measurements of metallicity, I as- 
signed [Z] — 0.0, with an RMS uncertainty of 0.2 
dex. These parameters fairly accurately character- 
ize the distribution of EBs having observational 
estimates of [Z], though the mean metallicity of 
stars in the solar neigh borhood is smaller tha n 



this by about 0.17 dex (jNordstrom et al.M2004f) . 



Effective temperatures for the TAG stars came al- 
most entirely from photometry; for stars cooler 
than 8000 K the estimated uncertainties in Tcs 
were typically between 50 K to 250 K, though for 
some hotter stars the uncertainties were as large as 
800 K. Data for the 156 stars used in this analysis 
are given in Table 1. 

For each star, I applied the MCMC procedure 
with the search over stellar models constrained by 
the star's mean density (computed from the TAG 
mass and radius), Tcs, and [Z]. The uncertainties 
I assumed in these quantities are described above. 
The MCMC/amoeba procedure then yielded new 
estimates of the stellar mass, radius, luminosity, 
age, mean density, Toff, and [Z], where the proce- 
dure necessarily returned values of the last 3 quan- 



tities that were much less than 1 standard 

deviation from those that were provided as in- 
put. The values of M*, i?*, and A^, derived by 
this method are listed in the last 3 columns (la- 
beled "fit") of Table 1, and in the corresponding 
columns of Table 2. 

To compare input and output masses and radii, 
for each star I computed the normalized discrep- 



ancy 5M^ 



{Moutput — M input) I M input, and 



i5i?* = {Routput - Rinput)/ Rinput- Figure 1 shows 
SRf, plotted against the stellar masses estimated 
by TAG. Since the fits included a strong constraint 
on the stellar mean density, (as calculated from 
the provided values of A'/* and i?*) I obtained the 
nearly perfect equality 5M^, = 2,6 R^, . One of the 
quantities (5M* or 5R^ is therefore redundant, and 
can easily be computed from the other. In what 
follows, I will frame the discussion in terms of the 
discrepancy in radius (5i?, . 

From the Figure, it is evident that the p.^, 
method of estimating stellar radii works fairly well. 
Aside from a few outliers, the RMS deviation of 
is 2.6%, while that of 5M^ (not plotted) is 
3 times as large at 7.8%. The discrepancies are, 
however, predominantly negative, and they vary 
systematically with stellar mass. The sense of this 
trend is that the masses and radii derived using 
the p* method are generally too small, and in per- 
centage terms, the more so for smaller stars. I 
have plotted discrepancies as a function of the stel- 
lar mass because this parameter is independent of 
stellar age. But since almost all of the stars on 
this plot are on the main-sequence or only slightly 
evolved, changing the independent variable to stel- 
lar radius or Teff would show essentially the same 
picture. These parameters are so tightly corre- 
lated along this part of the main sequence that 
without further information, one can hardly say 
which of them (or what combination of them) is 
most central to the observed effect. 

Fitting 5R^: to a linear function of log yields 

5R. = (-0.0175±0.0023)+ (0.0561±0.0076)logM* 

(1) 

This fit shows unambiguously that 5R^, is negative 
and has a trend with mass. The coefficients are 
nonzero with high statistical significance - more 
than 6cr for each, with the fits removing about one- 
third of the variance in the data set, and shrinking 
the statistic by 40. Of course, one can also 
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Fig. 1. — The relative radius discrepancy Ji?* for 

stars in the TAG sample, plotted against stellar 
mass . The error bars show 1-cr (assumed sym- 
metrical) errors derived from the MCMC analysis, 
added in quadrature to the radius errors quoted by 
TAG. The stars for which there are individual es- 
timates of [Z] are shown as filled circles. I show 
the Sun's position as a large symbol, the dot- 
ted line show the results of a fit to SR^ as a linear 
function of logM*. The dashed line is a fit to 
the picccwisc linear function of logM.^ described 
in the text. Stars excluded from these fits are plot- 
ted with special symbols, and are described in the 
text. 



fit more complicated functions to the data. For 
instance, a piecewise linear function that has zero 
slope for M, > IAMq fits the data marginally 
(but not compcUingly) worse than docs a single 
linear function; this is shown as dashed lines in 
the Figure. It is described by 

6R^ = (0.0980 ± 0.020) log(Af./1.4), M, < 1.4 

(2) 

= (-0.0021 ±.098), > 1.4 

I excluded 6 outlying stars from the least- 
squares fits in Fig. 1. These are the two bina- 
ries V1174 Ori and OGLE 051019, as well as the 
primary components of SZ Ccn and TZ For. All 
of these stars except perhaps SZ Gen are unusual 
as regards their evolutionary state or metallicity. 
Their individual peculiarities are described in sec- 
tion 4.4 below. 

So far I have discussed only the masses and 
radii of stars in the TAG sample. Considering 
the best-fitting age values gives further evidence 
that the evolution models are failing for at least 
some stars. Stars with masses larger than 1 Mq 
have, almost without exception, inferred ages that 
are less (usually much less) than that of the Milky 
Way. On the other hand, stars below solar mass 
have, with very few exceptions, inferred ages con- 
siderably greater than the Sun's, and most often 
greater than a Hubble time. The record-holder 
in this regard is VI 174 Ori B, with an apparent 
age that is greater than 100 GY, but many other 
low-mass stars display impossible ages. 

The likely reason for these excessive age es- 
timates is that the low-mass stars in the TAG 
sample have mean densities that arc so low that 
stars with the observed values of Tes and [Z] at- 
tain them only after they are significantly evolved. 
Given the extended main-sequence lifetimes of 
these stars, such evolution can take a very long 
time. A related possibility is that, for low-mass 
stars, the measured effective temperatures are sys- 
tematically too low. If this were the case, then the 
ZAMS model for a cool star would, in addition to 
being too cool, be smaller and more dense than it 
should be. Getting such a model to the observed 
p* would then require greater age than it should, 
and perhaps greater than possible. In either case, 
one can be sure that model-fitting solutions that 
require unphysical ages indicate something wrong 
with the models. 
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4. Discussion 

Systematic difFcicinccis between YY models and 
the samples of stars described in the previous sec- 
tion are concerning, even if they are rather small. 
But the way in which the differences are inter- 
preted (and perhaps corrected) depends on how 
they arise. Three general sources of systematic 
error seem possible. First, the masses and radii 
estimated by the method may be incorrect be- 
cause of systematic errors in the input data. Sec- 
ond, the chosen samples of stars may be systemat- 
ically different from the galactic mean, or at least 
from those stars that one expects to encounter in 
planet searches. Last, the evolution models may 
simply be erroneous in some respects. Any or all 
of these possibilities may be true to some extent. 
I consider each of them below. 

4.1. Systematic Errors in Input Data 

Biased input data may of course lead to biased 
estimates of stellar mass and radius. In particu- 
lar, values of Tes and [Z] derived from compar- 
isons between observed spectra and stellar atmo- 
spheres models may not correspond exactly to Toff 
and [Z] as used in the stellar evolution models. 
Moreover, such inaccurate correspondence might 
plausibly depend on Tgg itself, which in principle 
could generate the observed lack of agreement, (p* 
might be in error also, but as applied to the TAG 
sample of EBs and to the seismically-constrained 
stars it is correct by definition, and as applied to 
transiting planets it depends only on such well- 
established physics that its reliability seems as- 
sured.) 

Systematically inaccurate reddening estimates 
might also influence the masses and radii, by gen- 
erating systematic errors in the observed T^g val- 
ues. This effect would however be strongest for 
the hottest and most luminous stars in the sam- 
ple, since such stars tend to be the most reddened. 
There is no evidence in Fig. 1 for an effect of this 
sort, so I have not pursued this possibility. 

To assess whether biased Tcff or [Z] values 
might account for the observed discrepancy, I 
estimated the sensitivities dlogR^/dlogTes and 
d log R^,/d[Z] by perturbing the input values for 
the TAG sample. Figure 2 shows these sensitivi- 
ties as a function of logM*, computed using per- 
turbations of 0.007 in logTeff and 0.1 dex in [Z]. 




0.0 0.2 0.4 

log(M) (solar units) 



0.0 0.2 0.4 

log(M) (solar units) 



Fig. 2. Sensitivities of the estimate of i?* re- 
sulting from the p* method for stars in the TAG 
sample of eclipsing binaries. The sensitivity to [Z] 
(upper panel) was computed using 5[Z] = 0.1; the 
sensitivity to T^s (lower panel) used STes = 0.007. 
Dashed lines show the results of robust fits to lin- 
ear functions of log M* . 



a 
s 





Fig. 3. — p* vs logTeff plot for the stars in the 
TAG sample. Yonsei-Yale evolutionary tracks for 
stars with = {0.75, 1., 1.4, 2., 5.}Mq arc over- 
plotted as solid lines. Loci of constant central Y 
are plotted as dashed lines, corresponding to the 
ZAMS, (solar), and main-sequence turnoff. The 
Sun is indicated with a © symbol. 
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Not all of the binaries yielded valid differences, 
because of jumps between multiple solutions as de- 
scribed in section 2. For the majority of stars, 
however, the results indicate that the sensitivi- 
ties are only weakly mass-dependent. A robust 
straight-line fit to the sensitivities yields 

91ogi?,/91ogrcfr = 0.48 - 0.12 log (3) 

d\ogR^ld[Z] = 0.048 - 0.029 log . (4) 

Both of these sensitivities are well defined for 
logM* > 0.15, and less so for smaller masses. 
Even at 1 solar mass, however, it seems safe to esti- 
mate that to cause the observed -1.7% discrepancy 
in would require systematic errors of about 
200 K at TcH = 5800 K, or more than 0.15 dex in 
[Z\. Systematic errors of this size are essentially 
impossible for Toff, and implausible for [Z]. Thus, 
while one cannot altogether discount systematic 
errors in these observed parameters as contribut- 
ing to the discrepancy between TAG parameters 
and those from the p* analysis, it seems unlikely 
that they make an important contribution. 

4.2. Bias in the TAG Sample 

Perhaps the most striking point to be made 
from Fig. 1 is that the Sun, to which the YY 
models are calibrated, is not closely representa- 
tive of the stars of its general type in this sample. 
Compared to its neighbors in Fig. 1, it is assigned 
a larger mass and radius than is typical, by about 
5% and 1.6%, respectively. These discrepancies 
clearly do not arise from errors in the Sun's prop- 
erties. There seem to be three possibilities: the 
EB stars that fall near the Sun in mass and ra- 
dius are subject to systematic errors of a few per- 
cent in their mass and radius measurements, or 
they are subject to selection bias of similar mag- 
nitude, or there is genuine star-to-star variability 
in the stellar mass-radius-luminosity relationship. 
Again, none of these possibilities excludes the oth- 
ers; all may be acting to some degree. 

To address these issues it is helpful to display 
the data in a way that separates the observational 
data (EB masses and radii) from issues of evolu- 
tion modeling and data fitting. Figure 3 shows 
p* plotted against logTcff for the TAG sample of 
stars, and for the Sun. The axes are oriented as 
in a color-temperature diagram, with temperature 



increasing to the left, and (a proxy for lu- 
minosity) increasing upwards. Overlaid are (solid 
lines) YY evolution tracks for a sampling of stel- 
lar masses, and (dashed lines) loci correspond- 
ing to the zero-age main sequence (ZAMS) , main- 
sequence middle age (with about 71% helium in 
the core, corresponding to the Sun's evolution- 
ary state), and the age of turnoff from the main 
sequence. Even more than in Fig. 1, the Sun 
lies at or near one boundary of the region popu- 
lated by stars: it has as large or larger than all 
other stars with similar Tefj. Evolution in stars of 
roughly solar mass carries them upward on this di- 
agram, to lower p* and higher luminosity. So, bar- 
ring a statistical fluke, it seems either that all of 
the low-mass stars in the TAG sample are at least 
as evolved as the Sun, or that evolution models 
that accurately describes the Sun do not describe 
the average properties of the TAG sample. 

Two obvious sources of bias in the TAG sam- 
ple are Malmquist bias, which will favor detection 
of more luminous binaries, and a geometrical bias 
towards stars with larger radii, since (assuming 
random orbital inclinations) these are more likely 
to present themselves as EBs. These effects are 
however inadequate to account for the strength or 
the morphology of the observed bias. According 
to the YY evolution tracks, the present-day Sun is 
about 36% more luminous and 12% larger in ra- 
dius than it was on the ZAMS. Its detectability 
as the larger member of an EB should be propor- 
tional to the volume within which it is brighter 
than some limiting apparent magnitude, multi- 
plied by the ratio of solid angles within which the 
rotation axis must lie in order for eclipses to oc- 
cur. The first factor is at most proportional to 

3/2 

LJ , while for well-detached EBs such as found 
in the TAG sample, the second factor is propor- 
tional to ii*. Sun-like stars in the second half of 
their main-sequence lifetime thus should be over- 
represented among the TAG stars by a factor of 
at most about 1.8, relative to younger stars of 
similar mass. It is surprising to find the latter 
group totally absent. Moreover, the factors favor- 
ing detection of somewhat evolved stars should be 
similar for stars of all masses, whereas the trend 
seen in Fig. 3 is quite different. For stars with 
masses greater than 1.4 Mq, the ZAMS is coin- 
cident with the observed maximum-p* boundary, 
and for more massive stars there are roughly equal 
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numbers in the first and last halves of their main- 
sequence lifetimes. But for masses below 0.8 M©, 
no stars are found anywhere within the YY main 
sequence (although the significance of this obser- 
vation is shaky because of the small number of 
TAG stars in this mass range). On balance, it 
seems unlikely that Malmquist and geometrical se- 
lection bias by themselves provide an adequate ex- 
planation for the paucity of high-p.^ Sun- like stars 
in the TAG sample. 

Since most of the TAG stars lack metallicity 
determinations, another possible bias is that the 
true [Z] values for these stars differs systemati- 
cally from the assumed Gaussian distribution with 
[Z] = 0.0 ± 0.2. To test this idea I performed fits 
to the observed SR vs. logAf* dependence, us- 
ing only the 32 stars with measured metallicities. 
For both the linear and piecewise linear fits, the 
coefficients of the mass-dependent part of the vari- 
ation agreed with those in Eqs. (1-2) within Icr, 
and for both fits these coefficients differed from 
zero by about 2.5<t. Thus, the subset of stars with 
measured metallicities gives results that arc con- 
sistent with those from the full sample, but with 
worse precision. 

To explore a possible [Z] bias further, I artifi- 
cially altered the model [Z] values as a linear func- 
tion of log Teff, in order to remove the mass depen- 
dence seen in Fig 1. For all but the least massive 
(M < O.QAIq) stars, I achieved this by increasing 
[Z] by an amount S[Z] = +0.1-2 log(Teff/T0), but 
only for T^s < 6400K. This implies that to pro- 
duce the observed radii, the ~70 low-mass stars 
in the TAG sample would need to be systemati- 
cally metal-rich compared to the Sun; about 25% 
would require [Z] > 0.2, and about 10% would 
need [Z] > 0.4. While such a distribution of [Z] 
is conceivable, it is implausible a priori. A survey 
of about 12,000 bright stars in the solar neigh- 
borhood ( Nordstrdm et al] l2004[ ) shows a distri- 
bution that is centered at [Z] ~ —0.17, with a 
characteristic width of about 0.22 dex. Signifi- 
cantly metal-rich stars are therefore fairly rare; 
only 4% of Nordstrom's sample show [Z] > 0.2, 
and fewer than 0.4% have [Z] > 0.4. Since no 
obvious metallicity selection is operating in the 
choice of EBs in the TAG sample, it is implau- 
sible that the sample should be as skewed towards 
high [Z] as is required to explain Fig. 1. On the 
other hand, systematically high [Z] would help re- 



duce the large ages attributed to many low-mass 
stars, because increasing metallicity tends to de- 
crease stellar mean density at constant TeS and 
age. A final curiosity is that the very lowest-mass 
stars (smaller than about 0.6 Mq) show radius 
changes with [Z] that are smaller than (and some- 
times opposite in sign to) those of their more mas- 
sive brethren. Evidently it is dangerous to draw 
conclusions base on these stars, presumably be- 
cause there are many changes in important phys- 
ical processes as one nears the boundary between 
K and M stars. 

A bias that is more likely to be important re- 
lates to the stars' angula r momenta and mag- 
netic activity. Torres et al. (l2006f) . Lopez-Moraled 
(|2007l ). and [Morales et al.l (|2008[ ) have noted that 
some rapidly-rotating, magnetically-active EB 
component stars have temperatures that are lower 
and radii that are larger than expected for their 
masses or luminosities. Conspicuo us examples 



of su ch stars include V1061 Cyg (ITorres et al 



20061) . and CV Boo (iTorres et al.ll2009l). Both 
models (jChabrier et al.ll2007l : IClausen et al.ll2009f ) 

and general structural considerations suggest that 
magnetic processes in stellar outer convection 
zones are responsible for inflating stars in this 
way. Moreover, as pointed out by TAG and oth- 
ers, EB systems are strongly selected for small, 
short-period orbits. As a result, their components 
usually are tidally locked, and rotate rapidly com- 
pared to all but the youngest field stars. Thus, one 
expects the TAG sample to be biased in the sense 
that its members have larger radii and smaller 
Teff than would a similar sample of single stars. 
The paucity of high-density of stars seen in Fig. 
3 is consistent with this expectation. The lack of 
dense stars is most evident for sub-solar masses 
(for which stars have deep convection zones), and 
weak or absent for stars above 1.4 Mq, for which 
surface convection zones almost vanish. 

4.3. Asteroseismic Sample 

The 15 stars with properties listed in Table 2 
have measurements of their pulsation frequencies 
that constrain their large frequenc y separations 



Av, a nd hence the values of (e.g.. [Hansen et al 



(tiool ). When combined with other astronomical 
data (the details vary from star to star), this infor- 
mation allows accurate estimates of i?* and M*. 
In most cases, the stellar oscillations have been ob- 
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served in the radial velocity signal; for this to be 
feasible, the stars must have small values of the ro- 
tational line broadening vsini. Hence, these stars 
should not suffer from the rotation/activity bias 
ascribed to the TAG sample. Two stars, a Cen A 
and B, appear in both the TAG and the astero- 
seismic lists. 

Fig. 4 shows the result of applying the p* 
method to the asteroseismic sample of stars and 
plotting mass and radius discrepancies as in Fig. 
1. Expressed in the same terms as Eqns (1) and 
(2), a fit to these data yields 

SR^ = (0.0013±0.0068) - (0.002 ±0.046) log M* 

(5) 

and 

6R^ = (0.042 ± 0.075) log(M*/1.4), < 1.4 
SR^ ^ (0.005 ± .009), > 1.4 (6) 

Evidently the YY models fits the asteroseismol- 
ogy sample within the errors. Moreover, the zero 
points of the linear fits to the TAG and astero- 
seismology samples disagree by 2.5 a, so there is 
moderately strong evidence that the two distribu- 
tions are indeed drawn from different populations. 
The fits also yield discrepant slopes for the two 
samples, but the asteroseismology sample slope is 
so poorly constrained that the difference between 
slopes is only about 1 a. 

Figure 5 shows the vs Tes plot for the as- 
teroseismic sample. One star, ( Hyi, is a true gi- 
ant with logp* = —2.6, and hence falls outside 
the plotted range of the figure. There are 4 stars 
out of 15 in the asteroseismology group that have 
equal or greater density and equal or younger evo- 
lutionary age than the Sun. This is consistent with 
the hypothesis that low-mass EB component stars 
have systematically lower density than do slowly- 
rotating field stars of the same Teff. As yet, how- 
ever, the number of stars with asteroseismic mass 
and radius measurements is too small to justify a 
firm conclusion in this regard. Moreover, as noted 
earlier, asteroseismic parameter estimation has its 
own model dependencies. Agreement between as- 
teroseismic and /9* methods for estimating param- 
eters may say more about similarities between the 
models used than it says about the accuracy of 
either set of models in representing real stars. 
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4.4. Stars with Unusual Properties 

In Figures 1 and 4, a few stars stand out as be- 
ing significantly inconsistent with the mean rela- 
tion between SR^, and . I did not include these 
stars in the least-squares fits described in section 
3, so it is worthwhile to understand in what other 
ways they differ from the norm. In most cases 
there is a ready explanation for their odd radii. 

The two stars overplotted with asterisks in Fig. 
2 have fractional radius differences smaller than 
-0.1. They are the compon ents of VB1174 Ori . 
a pre-main-sequence system (IStassun et al 



20041) 



with oversize components that are still contract- 
ing. 

The two stars overplotted with tria ngles are 



comp onents of OGLE-05019.64-685812.3 (jPietrzvhski et 
20091 ). This member of the Large Magellanic 
Cloud is the faintest and the second most metal- 
deficient system in the TAG list of binaries (after 
V432 Aur). Its assigned metalhcity {[Z] = -0.5) 
has not been measured directly; rather it is in- 
ferred from the mean metallicity of its popula- 
tion. Its small radius relative to the YY models 
might therefore result if its metallicity were actu- 
ally about [Z] = —1.1, or as a result of errors in 
the models themselves, or in how the latter are 
used. The last explanation seems most likely, as 
these stars may in fact be clump giants, which lie 
beyond the oldest evolutionary stage tabulated in 
the YY models. 

The 2 stars with SR^, ~ -1-0.05 and overplotted 
with diamonds are the primary components of SZ 
Cen (with M, = 2.31) and TZ For (M* = 2.05). 
The SZ Cen system is poorly represented by evolu- 
tion trac ks for coeval stars wit h equal initial com- 
position (|Gronbech et al. I ll977l) . but the reason for 
this failure is obscure. TZ For consists of a pair of 
giants, with the pri mary component's ra dius be- 
ing more than 8 Rq ([Andersen et al ] ll99lh . These 
authors found that the more luminous, cooler com- 
ponent of this system (identified as TZ For A in 
the compilation by T AG, but as TZ For B in 
Andersen et"al] (|l99lf l) can be reconciled with the 
properties its companion only if the former is in 
the core helium burning stage of evolution. Since 
the YY models extend only to the tip of the giant 
branch, this star cannot be represented by them. 

Finally, the lowest-mass and most discrepant 
star in the list of stars with asteroseismic data is r 



Get ( Teixeira et al.|[2"009l) . This star is well known 
for having extremely low magnetic activity, yet its 
asteroseismic radius appears 5% larger than the 
YY models would indicate given its T^s and [Z], 
even if one allows the unreasonable age of 24 GY. 
In this case the asteroseismic data were compro- 
mised by instrumental noise, and hence may be 
open to question. Rather than speculate about 
causes, it seems best in this instance to wait for 
further observations. 

4.5. Conclusions 

The original motivation for this paper was to 
assess the accuracy of the p, method for estimat- 
ing radii of the host stars of transiting extraso- 
■^jjir planets. Since the reliability of the method 
IS almost entirely dependent upon the reliability 
of stellar evolution models, this question can also 
be framed in terms of the ability of such mod- 
els (in the present case, the Yonsei-Yale models) 
to represent the true mass-radius-luminosity de- 
pendencies of real stars. The short answer to the 
initial question is yes, without any corrections or 
adjustments the method implemented with the 
YY evolution tracks will yield stellar radii that 
are accurate within (at worst) about 5%, for a 
large majority of stars. This accuracy in radius is 
already good enough to support useful measure- 
ments of planetary density from transit measure- 
ments by Kepl er and other space borne photom- 
etry missions (jSeager et al. I l2007l) . although the 
corresponding 15% uncertainty in mass may be 
problematic. However, if one confines attention to 
the asteroseismic sample of single, slowly-rotating 
stars (Table 2), the accuracy of p^-based radius 
estimates appears considerably better, perhaps as 
good as 2% RMS. At this level, the accuracy of 
radius estimates is limited by probable system- 
atic errors not only in the evolution models, but 
to a similar degree by errors in estimates of stel- 
lar Teff and [Z]. Improving this situation will 
require further work to obtain accurate spectro- 
scopic estimates of these quantities and to under- 
stand the correlated errors between them, espe- 
cially for spectra with lines that are significantly 
broadened by stellar rotation. 

Although the method (with YY models) 
works well on average, it shows mass-dependent 
systematic errors of several percent when applied 
to components of eclipsing binary systems. The 
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sense of these errors is that observed EB compo- 
nent radii and masses are larger than those pre- 
dicted by the YY models; these differences are 
largest for low-mass stars. A plausible reading 
of the data (but one that is not compelled by 
statistics) is that the radius differences are zero 
for stars with more than about 1.4 solar masses, 
corresponding to the mass at which surface con- 
vection zones shrink to a negligible fraction of 
the stellar radius. For many years, it has been 
noted that rapidly-rotating, magnetically-active 
stars may have radii that are up to 10% large r than 
their masses and luminosities w ould imply (jHoxie 
19731 : lLacvlll977l : |Poppej[l997h . These claims are 
based for the most part on results for unequal- 
mass EB components, for which it proves impossi- 
ble to fit masses, radii, and luminosities if the com- 
ponents are assumed to be coeval and to have the 
same compos ition. Examples of such stars include 
CV Boo B iT orres et_alj l2009|) a nd V1061 Cyg 
B (iTorres et a l. 2006.). Theory (jChabrier et al 



20071 ) suggests that the radius expansion results 
from blockage of the convective energy flux by 
strong magnetic fields in the stellar convection 
zones. The mass dependence of EB radius dis- 
crepancies seen in Fig. 1 is thus likely to be the 
population-averaged expression of the same phys- 
ical process. 

It is poss ible that th e stellar models calculated 
by Baraffe (|Chabrier fc Baraffe 1997) fit the ob- 
served properties of cool, low-mass stars better 
than the YY models do. These models are in- 
tended to treat such stars, using non-gray at- 
mospheres and smaller mixing-length parameters 
that act to reduce the efficiency of energy trans- 
port in the outer stellar layers. It would be worth- 
while (but is outside the scope of this paper) to 
perform a detailed comparison of the TAG masses 
and radii with those inferred using the method 
and the Baraffe models. 

Stars with asteroseismic estimates of do not 
show the radius discrepancies found in the TAG 
sample. This sample of stars are all magnetically 
inactive, slow rotators. The good agreement be- 
tween their observed radii and those predicted by 
the YY models argues in favor of magnetic activity 
having a visible influence on the radii of low-mass 
stars. The number of asteroseismically-measured 
stars is, however, as yet too small to allow a firm 
conclusion that these stars are drawn from a differ- 



ent population than the eclipsing binaries. Future 
work will surely improve this situation, especially 
insofar as spaceborne photometry allows astero- 
seismic measurements on rapidly-rotating stars, 
which are difficult to measure with Doppler meth- 
ods. 

I am grateful to G. Torres for his valuable com- 
ments on an early draft of this paper, and to him, 
J. Andersen, and A. Gimencz for providing me the 
information on which Table 1 is based in advance 
of publication. I also thank an anonymous referee, 
whose suggestions materially improved the discus- 
sion in 64.2. 
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Table 1 

Eclipsing Binary Stars with Accurate Mass and Radius Estimates. 
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FL Lyr B 


0.958 





.012 


0.962 


0, 


.028 


5300 


100 


0.000 


0. 


.200 


0.815 


0.913 


15.25 


V432 Aur A 


1.204 





.006 


2.430 


0. 


,023 


6080 


85 


-0.600 


0, 


.050 


1.168 


2.405 


4.36 


V432 Aur B 


1.079 


0, 


.005 


1.224 


0, 


.007 


6685 


85 


-0.600 


0, 


.050 


1.071 


1.220 


3.55 


EW Ori A 


1.174 


0, 


.012 


1.134 


0. 


Oil 


5970 


100 


0.000 


0, 


.200 


1.074 


1.100 


3.67 


EW Ori B 


1.124 





.009 


1.083 


0, 


Oil 


5780 


100 


0.000 


0. 


.200 


0.989 


1.038 


6.03 


BH Vir A 


1.166 





.008 


1.247 


0, 


.024 


6100 


100 


0.000 


0. 


.200 


1.162 


1.246 


3.12 


BH Vir B 


1.052 





.006 


1.135 


0, 


,023 


5500 


200 


0,000 


0. 


.200 


0.929 


1.086 


11.04 


ZZ UMa A 


1.139 





.005 


1.514 


0, 


,019 


5960 


70 


0,000 


0. 


.200 


1.148 


1.516 


5.54 


ZZ UMa B 


0.969 


0, 


.005 


1.156 


0. 


.010 


5270 


90 


0.000 


0, 


.200 


0.846 


1.103 


18.47 


a Cen A 


1.105 


0, 


.007 


1.224 


0, 


,003 


5824 


26 


0.240 


0, 


.040 


1.124 


1.231 


4.69 



Table 1 — Continued 



Star 


M, 


M, 


i?„ 


R, 


Toff 


Tcff 




[. 




M 


fit 




Agojit 




(Mo) 


RMS 


(Rq) 


RMS 


(K) 


RMS 






RMS 


(Mo) 


(Rq) 


(GY) 




0. 


934 





.006 





.863 





.005 


5223 


62 


0, 


,240 





,040 


0, 


.926 


0.859 


3 .79 


iN^jl^loo IVrt VIZ J\ 


1, 


103 





.007 


1 


.426 





.019 


5900 


100 


-0, 


100 


0, 


090 


1, 


.043 


1.401 


1 .O 1 


iN^jL^ioo J\.rt Viz -tS 


1, 


081 





.007 


1 


.374 





.019 


5875 


100 


-0, 


100 


0, 


090 


1, 


.052 


1.358 


7.03 


V oDo i^yr A 


1, 


074 





.008 


1 


.400 





.016 


5665 


100 


0, 


,400 


0, 


100 


1, 


.200 


1 .452 


4,91 


V ODo i^yr 15 


0. 


827 





.004 





.768 





.006 


4900 


100 


0, 


,400 





,100 


0. 


.861 


.779 


3.04 


V uoD oen rt 


1, 


052 





.005 


1 


.019 





.004 


5900 


85 


-0, 


,200 





,080 


0. 


.958 


.988 


5.42 


V ooD ocn JiS 


0. 


854 





.003 





.830 





.004 


5000 


100 


-0, 


,200 





,080 


0. 


.681 


0. 771 


25.55 


V JjOO J\ 


1, 


032 





.013 


1 


.263 





.023 


5760 


150 


0, 


,000 





,200 


1, 


.029 


1.262 


7.91 


V l500 D 


0. 


968 





.012 


1 


. 174 





.023 


5670 


150 


0, 


,000 





,200 


0. 


.925 


1 .160 


1 1,80 


1 TA r\-r-\ a 
V i i / 4 kjri rt 


1, 


006 





.013 


1 


.338 





,011 


44 ( U 


120 


0, 


,000 


0, 


200 


0. 


.698 


1. 184 


42,39 


■\7'1 1 7/1 R 
V i i ( 41: KJTl ±5 


0, 


,727 


0, 


.010 


1 


.063 


0, 


,011 


3615 


100 


0, 


,000 


0, 


,200 


0, 


.428 


0.891 


108,00 




0, 


,983 


0, 


.008 


1 


.110 


0, 


,023 


5780 


100 


0, 


,000 


0, 


,200 


0, 


.981 


1 .108 


7 7Q 


U V r^SC 15 


0, 


,764 


0, 


.004 





.835 


0, 


,018 


4750 


SO 


0, 


,000 


0, 


,200 


0, 


.710 


0.813 


28,42 


CG Cyg A 


0, 


,941 


0, 


.014 





.893 


0, 


,012 


5260 


180 


0, 


,000 


0, 


,200 


0, 


.834 


0.858 


11.18 


CG Cyg B 


0, 


,815 


0, 


.013 





.838 


0, 


,011 


4720 


60 


0, 


,000 


0, 


,200 


0, 


.700 


0.795 


29.00 


RW Lac A 


0, 


,926 


0, 


.006 


1 


.187 


0, 


,004 


5760 


100 


0, 


,000 


0, 


,200 


0, 


.973 


1.210 


9.63 


RW Lac B 


0, 


,869 


0, 


.004 





.964 


0, 


,004 


5560 


150 


0, 


,000 


0, 


,200 


0, 


.905 


0.977 


9.64 


HS Aur A 


0, 


,898 


0, 


.019 


1 


.005 


0, 


,024 


5350 


75 


0, 


,000 


0, 


,200 


0, 


.834 


0.980 


15.82 


HS Aur B 


0, 


,877 


0, 


.017 





.874 


0, 


,024 


5200 


75 


0, 


,000 


0, 


,200 


0, 


.813 


0.852 


13.09 


GU Boo A 


0, 


,609 


0, 


.006 





.627 


0, 


,016 


3920 


130 


0, 


,000 


0, 


,200 


0, 


.514 


0.592 


56.22 


GU Boo B 


0, 


598 


0, 


.006 





.623 


0, 


,016 


3810 


130 


0, 


,000 


0, 


,200 


0, 


.493 


0.584 


61.11 


YY Gem A 


0, 


,599 


0, 


.005 





.619 


0, 


,006 


3820 


100 


0, 


,000 


0, 


,200 


0, 


.498 


0.581 


57.80 


YY Gem B 


0, 


,599 


0, 


.005 





.619 


0, 


,006 


3820 


100 


0, 


,000 


0, 


,200 


0, 


.498 


0.581 


57.81 


CU Cnc A 


0, 


,435 


0, 


.001 





.432 


0, 


,005 


3160 


150 


0, 


,000 


0, 


,200 


0, 


.425 


0.434 


1.40 



^Data in all columns except the last 3 and in some cases [Z] and [Z]rms (as detailed below) are from pTbrres et al.| J2009D (TAG) 
^For stars with no observational estimate of [Z], I have set [Z] — 0., and [Z]ft.MS — 0-2, 



Table 2 

Stars with M* and i?* estimates based on asteroseismology. 



star 


M. 


M, 








Totf 


[Z] 


[Z] 


Mfit 


RfH 


Ago Jit 


Source 






(Af©) 


RMS 


(Rq) 


RMS 


(K) 


RMS 






RMS 




(-Ro) 


(GY) 






HD49933 


1 


.079 


0, 


.073 


1.385 


0. 


,031 


6650 


75 


-0 


,440 


0. 


,060 


1.150 


1.414 


3.47 


^^^^^ 1 


'2009) 


HD175726 


0, 


.993 


0, 


.060 


1.014 


0. 


,035 


6060 


50 


-0, 


,100 


0. 


,060 


1.063 


1.037 


2.34 


Bruntt 1 


^2609) 


JrliJlol'i^U 


1 


.311 


0, 


.063 


l.OifO 


0. 


,032 


DDZU 


1 nn 

lUU 


0, 


,000 


0. 


120 


1 "^Ql 

1 .oyi 




i.yo 


rsriniiT] \ 


'?nnQ^ 
^uuyj 


HD181906 


1 


.144 


0, 


.119 


1.392 


0. 


,054 


6365 


120 


-0, 


.110 


0. 


140 


1.224 


1.423 


3.10 


3^^^^ 1 


^2009) 


70 Oph A 


0, 


.895 


0, 


.005 


0.863 


0. 


,002 


5322 


40 


0, 


.011 


0. 


,050 


0.870 


0.855 


7.56 


^^^^^t al. 1 


:2008J 


a Ccn A 


1 


.105 


0, 


.007 


1.224 


0. 


,003 


5824 


26 


0, 


.240 


0. 


,040 


1.124 


1.231 


4.69 


^o^^^ct al. 1 


:2009) 


a Ccn B 


0, 


.934 


0, 


.006 


0.863 


0. 


,005 


5223 


62 


0, 


.240 


0. 


,040 


0.926 


0.859 


3.79 


^OTrcs ct al. 1 


:20095 


/i Ara 


1 


.100 


0, 


.010 


1.353 


0. 


,010 


5800 


90 


0, 


.300 


0. 


,050 


1.190 


1.389 


4.54 


Soriano & Vauclair 1 


:2009) 


Procyon 


1 


.497 


0, 


.037 


2.067 


0. 


,028 


6530 


90 


-0, 


.050 


0. 


,030 


1.477 


2.058 


2.24 


Eeeenbergcr ct al. 1 


'2005) 


rj Boo 


1 


.700 


0, 


.050 


2.790 


0. 


,040 


6030 


90 


0, 


.360 


0. 


,050 


1.715 


2.798 


1.96 


Garrier et al. 1 


'2005) 


/3 Hyi 


1 


.070 


0, 


.030 


1.814 


0. 


,017 


5872 


44 


-0, 


.030 


0. 


,050 


1.168 


1.867 


5.82 


North et al. 1 


'2007) 


S Eri 


1 


.215 


0, 


.020 


2.328 


0. 


,010 


5074 


60 


0, 


.130 


0. 


,080 


1.234 


2.344 


5.61 


Thevenin et al. ( 


f2005) 


CHyi 


2, 


.650 


0, 


.020 


10.300 


0. 


,100 


5010 


100 


-0, 


.040 


0. 


120 


2.777 


10.462 


0.47 


Thevenin et al. ( 


:2005) 


/3 Vir 


1 


.413 


0, 


.061 


1.703 


0. 


,022 


6059 


49 


0, 


.140 


0. 


,050 


1.301 


1.656 


3.65 


North ct al. 1 


^200^ 


T Get 


0, 


.783 


0, 


.012 


0.793 


0, 


,004 


5265 


100 


-0, 


.500 


0. 


,030 


0.660 


0.751 


24.29 


Tcixcira ct al. 1 


:2009) 



